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Abstract 

The  complexity  of  the  solid  oxide  fuel  cell  requires  the  use  of  very  large  computers.  For  this  reason,  models  based  on  various  engineering 
assumptions  have  been  created  so  that  the  numerical  calculations  may  be  performed  on  small  personal  computers. 

The  Hammerstein  models  are  special  kinds  of  nonlinear  systems  where  the  nonlinear  block  is  static  and  is  followed  by  a  linear  system. 
This  method  is  applied  to  model  the  static  and  dynamic  characteristics  of  the  solid  oxide  fuel  cell. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  solid  oxide  fuel  cell  (SOFC)  is  viable  for  generat¬ 
ing  electricity  from  hydrocarbon  fuels.  The  high  operating 
temperature  from  600-1000  °C  allows  internal  reforming  and 
promotes  rapid  kinetics  with  nonprecious  materials.  The  high 
temperature  of  the  SOFC,  however,  places  stringent  require¬ 
ments  on  its  materials. 

Padulles  et  al.  [1]  develop  a  SOFC  model  which  includes 
species  dynamics,  but  it  does  not  consider  temperature 
dynamics.  Hall  and  Colclaser  [2]  modeled  a  3-kW  SOFC 
but  they  did  not  take  into  account  dynamics  of  the  chemi¬ 
cal  species.  Achenbach  develops  a  mathematical  model  of  a 
planar  SOFC,  which  concentrates  on  effects  of  temperature 
changes  on  output  voltage  response  [3].  Temperature  dynam¬ 
ics  is  modeled  in  a  three-dimensional  (3D)  vector  space.  The 
same  author  investigated  the  transient  behavior  of  a  stand¬ 
alone  SOFC  caused  by  a  load  change  in  Ref.  [4].  It  shows 
that  the  relaxation  time  of  the  output  voltage  is  highly  related 
to  the  effect  of  temperature  dynamics. 
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The  application  of  the  autoregression  with  exogenous  sig¬ 
nal  (ARX)  identification  algorithm  to  compute  linear  system 
models  is  presented  in  Ref.  [5].  In  this  paper,  the  Hammer¬ 
stein  nonlinear  system  approach  is  used  for  identification  of 
the  nonlinear  system  model  of  SOFC. 

The  Hammerstein  models  are  special  kinds  of  nonlinear 
systems  where  the  nonlinear  block  is  static  and  is  followed 
by  a  linear  system.  These  models  have  applications  in  many 
engineering  problems,  and  therefore,  identification  of  Ham¬ 
merstein  models  has  been  an  active  research  area  for  many 
years.  There  exist  a  large  number  of  research  papers  in  the 
literature  on  the  topics  of  Hammerstein  model  identifications. 

Existing  methods  for  Hammerstein  model  identification 
can  be  roughly  divided  into  four  categories:  the  iterative 
method  [6,7],  the  over-parametrization  method  [8-11],  the 
stochastic  method  [12-14]  and  the  separable  least- squares 
method  [15].  In  most  cases,  the  structure  of  the  nonlinearity 
is  assumed  to  be  known.  Otherwise,  identification  becomes 
a  structural  estimation  problem  [16]. 

Due  also  to  the  particular  parameterization  used  in  Ref. 
[8],  an  extension  of  the  results  to  the  multiple  input,  multiple 
output  (MIMO)  setting  does  not  seem  to  be  straightforward. 
Gomez  and  Baeyens  [17]  proposed  a  noniterative  algorithm 
for  the  identification  of  Hammerstein  models,  which,  in 
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Nomenclature 

q 

heat  capacity  of  the  cell  unit  (J  (kg  K)-1) 

E0 

standard  reversible  cell  potential  (V) 

F 

Faraday’s  constant  (C  kmol-1) 

Ki  .  K " 

anode  (cathode)  inlet  partial  molar  enthalpies 
(J  mol-1) 

M 

partial  molar  enthalpies  at  stack  temperature 
(J  mol-1) 

i 

stack  current  density  (Am-2) 

A 

limiting  current  density  (Am-2) 

A 

exchange  current  density  (Am-2) 

/ 

stack  current  (A) 

Ki 

valve  molar  constants 

Kr 

constant  dependent  on  Faraday’s  constant 
and  number  of  electrons  in  the  reaction 
(kmol  (s  A)-1) 

M8 

mass  of  the  cell  unit  (kg) 

n 

number  of  electrons  participating  in  the  reac¬ 
tion 

Nh2 

hydrogen  flow  that  reacts  (kmol  s-1) 

N?,  Nf 

molar  flow  rates  (mols-1)  of  the  ith  reactant 
at  the  cell  input  and  output,  respectively 

Nf 

reaction  rate  (mol  s-1)  of  the  ith  reactant 

No 

number  of  cells  in  stack 

Nm  Nm 

anode  (cathode)  total  inlet  molar  flow 
(mol  s-1) 

P 

cell  pressure  (atm) 

P  dc 

stack  dc  power  (W) 

r 

ohmic  resistance  (Q  m2) 

R 

gas  constant  (8.31  J  (mol  K)-1) 

P'dh  Pci 

anode  (cathode)  total  rate  of  production  of 
species  (mol  s-1) 

r 

stack  solid  average  temperature  (K) 

7° 

temperature  constant  (K) 

v 

compartment  volume  (m3) 

Vdc 

cell  voltage  (V) 

Ve 

volume  of  the  cell  unit  (m3) 

V0 

open-circuit  reversible  potential  (V) 

rin  in 
Aa i  ’  Ac i 

anode  (cathode)  inlet  mole  fractions 

mole  fractions  of  species 

Greek  letters 

a 

electron  transfer  coefficient  of  the  reaction  at 

the  electrode 

O'r,  fir 

ohmic  resistance  constants 

7?  act 

activation  losses  (V) 

Ocon 

concentration  losses  (V) 

*h2 

time  constant  associated  with  the  hydrogen 
flow  and  is  a  function  of  temperature  (s) 

total  gas  components  in  anode  or  cathode 

contrast  to  [8],  applies  also  to  multivariable  systems,  allows 
a  more  general  representation  (using  base  functions)  for  the 
static  nonlinearity,  and  where  the  consistency  of  the  estimates 
is  guaranteed  even  in  the  presence  of  colored  output  noise. 
As  in  Ref.  [8],  the  main  computational  tools  employed  by  the 
algorithm  are  least- squares  estimation  (LSE)  and  singular 
value  decomposition  (SVD),  which  results  in  numerical 
robustness  under  weak  assumptions  on  the  persistency  of 
excitation  of  the  inputs.  Key  on  the  derivations  of  the  results 
in  Ref.  [17]  is  the  use  of  base  functions  for  the  representation 
of  the  linear  and  the  nonlinear  blocks  in  the  Hammerstein 
model. 

The  paper  is  organized  as  follows.  In  Section  2,  general 
principles  of  SOFC  are  explained.  The  multivariable  Ham¬ 
merstein  model  is  introduced  and  the  identification  problem 
is  formulated  in  Section  3.  Simulation  examples  illustrat¬ 
ing  the  performance  of  the  algorithms  on  SOFCs  are  pre¬ 
sented  in  Section  4,  and  finally,  conclusions  are  provided  in 
Section  5. 


2.  Solid  oxide  fuel  cell  dynamic  model 

The  proposed  stack  model  is  based  on  the  following 

assumptions: 

(1)  Stack  is  fed  with  hydrogen  and  air,  therefore,  the  fuel 
processor  dynamics  is  not  included. 

(2)  A  uniform  gas  distribution  among  cells  is  assumed,  since 
there  is  a  small  deviation  of  the  gas  distribution  among 
the  cells. 

(3)  There  is  no  heat  transfer  among  cells.  Each  cell  has  the 
same  temperature  and  current  density  [1]. 

(4)  The  channels  that  transport  gases  along  the  electrodes 
have  a  fixed  volume,  but  their  lengths  are  small,  so  that 
it  is  only  necessary  to  define  one  single  pressure  value  in 
their  interior.  The  ratio  of  pressures  between  the  interior 
and  exterior  of  the  channel  is  large  enough  to  consider 
that  orifice  is  choked  [1,18]. 

Fig.  1  shows  the  dynamic  model  of  SOFC,  along  with  its 

major  chemical  reactions. 

2.7.  Electrochemical  model 


The  change  in  concentration  of  each  specie  that  appears 
in  the  chemical  reactions  can  be  written  in  terms  of  input 
and  output  flow  rates  and  exit  molarity  due  to  the  following 
chemical  reaction  [19,20]: 


V  d 

Pi  =  N f  -  Nf  -  Nl 


RTdt 


In  agreement  with  the  basic  electrochemical  relationships, 
the  molar  flow  that  reacts  can  be  calculated  as: 


N\  = 


N0 1 


2KrI 


2  F 


(2) 
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Fig.  1.  Dynamic  model  of  SOFC. 


The  cell  utilization  is  defined  through  the  input  and  output  The  energy  balance  equation  for  each  cell  unit  is  as  fol- 

hydrogen  flow  rates  as  follows:  lows: 


u  = 


*i?2  ~  N°H2 


For  orifice  that  is  choked  [18],  molar  flow  of  any  gas  through 
the  valve  is  proportional  to  its  partial  pressure  inside  the  chan¬ 
nel  according  to  the  following  expressions  [1]: 


Nh2 

m2 


Nh2o 

Pn2o 


=  Ku2o 


Considering  the  hydrogen  partial  pressure, 


MSCS 


d  r 


p  dr 


=  N , 


m 

ta 


E 

i=  1 


m  / 1  in  i  s 


K) 


i=  1 


in 


+Ntc 


E 

i= 1 


in  zr^in  r  s 


K) 


y  ^  hi  R&  Pdc 

(8) 


i=i 


Under  the  ideal  gas  supposition,  the  partial  molar  enthalpies 
are  calculated  using 


^jtPK2  =  <-NH2-2KrI  (5) 

Applying  the  Laplace  transformation  to  the  above  equations 
and  isolating  the  hydrogen  partial  pressure,  yields  the  follow¬ 
ing  expressions: 


**  =  rrt  (<  - 2KJ) 

V 

th9  = - 

2  Kh2RT 


(6) 

(7) 


hi  =  hf  +  [  Cpj(u)  d u  (9) 

'l  Jr  of 

and  coefficients  of  the  specific  heats 

cp,/  —  +  hiT  +  c/T2  +  d[T 3  (10) 

are  encountered  in  standard  reference  tables.  The  reference 
enthalpy  stands  for  energy  at  standard  reference  temperature 
and  considers  the  heat  of  formation  for  each  gas  species  to 
account  for  energy  change  on  chemical  reaction. 


2.2.  Thermal  model 


2.3.  Nernst’s  equation 


The  fuel  cell  power  output  is  closely  related  to  the  tem¬ 
perature  of  the  cell  unit.  The  heat  storage  in  the  thin  fuel  unit 
gas  or  oxidant  gas  layer  is  neglected.  The  thin  fuel  unit  or 
oxidant  gas  layers  are  lumped  to  the  cell  unit  and  gas  layers 
are  assumed  to  have  the  same  temperature  as  the  cell  unit 
[2-4]. 


Applying  Nernst’s  equation  and  Ohm’s  law  (taking  into 
account  ohmic,  concentration,  and  activation  losses),  the 
stack  is  connected  in  series  and  the  stack  output  voltage  is 
represented  as  follows  [19-21]: 

Ldc  —  L0  rl  Tj&ci  7?con 


(ID 
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/  RT 

v0  =  n0Ie0  +  — 


In 


fHzfoJ 

*h2o 


(12) 


For  the  reason  that  the  reactant  is  consumed  at  the  elec¬ 
trode  by  electrochemical  reaction,  there  is  a  loss  of  potential 
due  to  the  inability  of  the  surrounding  material  to  maintain 
the  initial  concentration  of  the  bulk  fluid.  Concentration  loss 
equation  is  given  by  [2,19-21]: 


RT  (  i\ 

7? con  —  ~  In  I  1  —  T“  ) 

nF  \  iLJ 


(13) 


Activation  polarization  is  existent  when  the  rate  of  an  elec¬ 
trochemical  reaction  at  an  electrode  surface  is  controlled  by 
sluggish  electrode  kinetics  [20,21].  Activation  loss  equation 
is  as  follows: 


RT  (  i  \ 

^act=—  In  -  (14) 

ant  \i  o J 

a  is  the  transfer  coefficient,  which  is  considered  to  be  the 
fraction  of  the  change  in  polarization  that  leads  to  a  change 
in  the  reaction  rate  constant  and  its  value  is  usually  0.5  for 
the  fuel  cell  application. 

Tafel  plots  provide  a  visual  understanding  of  the  activa¬ 
tion  polarization  of  a  fuel  cell.  They  are  used  to  measure  the 
exchange  current  density,  given  by  the  extrapolated  intercept 
at  77 act  =  0  which  is  a  measure  of  the  maximum  current  that 
can  be  extracted  at  negligible  polarization,  and  the  transfer 
coefficient  (from  the  slope). 

The  usual  form  of  the  Tafel  equation  that  can  be  easily 
expressed  by  a  Tafel  plot  is 


77  act  =  a  +  b  In  i 


(15) 


where  a  =  (— RT/aNF )  In  z'o  and  b  =  RT/aNF. 

The  term  b  is  called  the  Tafel  slope  and  is  obtained  from 
the  slope  of  a  plot  of  77act  as  a  function  of  In  i.  The  Tafel  slope 
for  an  electrochemical  reaction  is  about  llOmV/decade  at 
room  temperature. 

Ohmic  losses  occur  because  of  resistance  to  the  flow  of 
ions  in  the  electrolyte  and  resistance  to  the  flow  of  electrons 
through  the  electrode  materials.  This  resistance  is  dependent 
on  the  cell  temperature  and  is  obtained  by  Ref.  [21] 


r  =  ar  exp 


(16) 


3.  Hammerstein  model  identification 

3.1.  Problem  formulation 

A  multivariable  Hammerstein  model  is  symbolized  in 
Fig.  2.  The  model  consists  of  a  zero-memory  nonlinear  ele¬ 
ment  Af(.)  in  cascade  with  a  linear  time  invariant  (LTI)  system 
with  transfer  function  matrix  G(q )  e  Hfxn( T).  It  is  accepted 
that  the  measured  output  yk  contains  an  unknown  additive 
noise  component  Vk> 


Fig.  2.  Multivariable  Hammerstein  model. 

The  input-output  relationship  is  then  given  by 
yk  =  G(q)Muk)  +  vk  (17) 

where  yk  £  Mm,  Uk  E  W1,  and  Vk  E  are  the  system  output, 
input,  and  measurement  noise  vectors  at  time  k,  respectively. 
It  will  be  assumed  that  the  nonlinear  block  can  be  specified 
as 

r 

Muk) =yfaigi(Mk)  (is) 

i—\ 

where  gf.)  =  W1  Mn(/  =  1,  . . . ,  r)  are  known  nonlin¬ 
ear  base  functions  and  at  e  Wlxn(i  =  1,  . . . ,  r)  are  unknown 
matrix  parameters.  Typically,  the  nonlinear  functions  gf.) 
are  polynomials  that  allow  the  representation  of  smooth  non- 
linearities,  but  they  can  also  be  radial  base  functions  (RBF) 
or  base  functions  generated  by  translations  and  dilations  of 
a  mother  function  (e.g.,  wavelets).  Any  smooth  function  in 
an  interval  can  be  represented  with  arbitrary  accuracy  by  a 
polynomial  of  sufficiently  high  order. 

From  another  point  of  view,  the  LTI  system  will  be  repre¬ 
sented  using  rational  orthonormal  bases  as  follows: 

p  i 

G{q)  =  YJblB,{q)  (19) 

1=0 

where  bieWnxn  are  unknown  matrix  parameters  and 
{J3i(q)}ff  0  are  rational  orthonormal  bases. 

The  identification  problem  is  to  estimate  the  unknown 
parameter  matrices  <zz-  eM”x”(/=  1 , ...  ,r)  and  bi  e  Wnxn 
(/  =  0,  1)  characterizing  the  nonlinear  and  the  lin¬ 

ear  parts,  respectively,  from  an  A-point  data  set  {up,  yk)k= i 
of  observed  input-output  measurements. 

3.2.  Identification  algorithm 

Substituting  Eqs.  (18)  and  (19)  in  (17),  the  input-output 
relationship  can  be  written  as 

p— l  r 

yk  =  'f2'2l  biajBi(q)gi(uk )  +  vk  (20) 

1=0  i= l 

From  Eq.  (20),  the  parametrization  (18)  and  (19)  is  not 
unique.  A  technique  that  can  be  used  to  obtain  uniqueness 
is  to  normalize  the  parameter  matrices  at  or  bn  for  example, 
assuming  that  |  \ai |  I2  =  1  or  |  \bi\  I2  =  1 .  Under  this  assumption, 
the  parametrization  is  unique. 

Specifying  now 

A  rP 

0  —  \b{^a  1 ,  . . . ,  b^Oy  •)  ...  1  b p —  ]  a  \  ^  . . . ,  b p —  ^  ay ]  (21) 
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Bo  (q)  gj(uk), Bp- 1  (q)  gj(uk ) 
Bp-\  ( q)gj(uk ), . . . , 


(22) 


Eq.  (20)  can  be  written  as 

yk  =  (P  <pk  +  vk  (23) 


where  0ap  has  been  defined  as 


alD0 

r-p  i 

a[  b 

a2  °0 

r-p  i 

a^b 

ajbj}  ajbj 


aTbT 

a\°P-\ 

aThT 

a2 °p- 1 


nThT 
ar  °p- 1 


(30) 


which  is  in  linear  regression  form.  Regarding  the  A-point  data 
set,  the  last  equation,  and  defining 


r  ■  i 

Yn  =  [yi,y2,  ■  ■  ■ ,  yn] 

(24) 

VN  =  [U1,  V2,  .  .  •  ,  VNf 

(25) 

&N  —  [01,  02,  •  •  •  ,  <Pn2 

(26) 

the  next  equation  can  be  written 

II 

^  H 

+ 

(27) 

A 

The  estimate  0  that  minimizes  a  quadratic  criterion  on  the 
prediction  errors  sN  =  Yjy  —  <PTN6  is  determined  by 

0  =  PnYn  (28) 

The  problem  is  how  to  calculate  the  parameter  matrices  <2; 

/V 

(i  =  1,  . . .,  r)  and  Z?/  (/  =  0,  1)  from  the  estimate  0  in 

(28).  From 

0  =  blockvec((9^)  (29) 

where  blocvec((9^)  is  the  block  column  matrix  obtained  by 
stacking  the  block  columns  of  0ap  on  top  of  each  other,  and 


where  0ap  =  abT  and  the  following  definitions  for  the  matri¬ 
ces  a  and  b , 

a  =  [a\,  ci2,  . . . ,  ar]T  (31) 

bHbl,b],...,bTp_1]T  (32) 

/V 

An  estimate  0 ^  of  the  matrix  0 ^  can  then  be  obtained 

✓V 

from  the  estimate  6  in  (28).  The  question  now  is  how  to  esti- 

/V 

mate  the  parameter  matrices  a  and  b  from  the  estimate  0ab-  It 

/V 

is  clear  that  the  closest  estimates  a  and  b  are  those  that  solve 
the  following  optimization  problem 

(a,  b)  =  arg  mini \\<S>a,b  -  (33) 

a,b  ^  ' 

The  solution  to  this  optimization  problem  is  supplied  by  the 

✓V 

SVD  of  the  matrix  0&p. 

The  structure  (17)— (19)  of  the  identified  Hammerstein 
model  can  be  interpreted  as  an  equivalent  LTI  model  whose 
inputs  are  the  actual  inputs  Uk  filtered  by  the  nonlinear  base 
functions  g/(.)  used  to  represent  the  static  nonlinearity  in  the 
Hammerstein  model.  This  interpretation  of  the  Hammerstein 
model  is  represented  in  the  block-diagram  of  Fig.  3. 


uk 


¥ 


¥ 


Uh 


Fig.  3.  Equivalent  LTI  model. 
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Table  1 


Operating  point  data 


Power  (kW) 

100 

Stack  voltage  (V) 

286.3 

Stack  current  (A) 

300 

Number  of  cells 

384 

Number  of  stacks 

1 

Open  circuit  voltage  for  each  cell  (V) 

0.935 

Input  fuel  flow  (kmol  s-1) 

1.2  x  icr3 

Input  air  flow  (kmol  s-1) 

2.4  x  10“3 

Cell  area  (cm2) 

1000 

Cell  temperature  (°C) 

1000 

Transfer  coefficient  a 

0.5 

Ohmic  resistance  constant  fiY 

-2870 

Ohmic  resistance  constant  aY 

0.2 

Temperature  constant  T°  (K) 

923 

Limiting  current  (Am-2) 

0.8 

ATht  (kmol  (atm  s)-1) 

8.43  X  10“4 

Fig.  4.  Output  stack  voltage  response  due  to  stack  current  step  increase. 


4.  Results 


The  identification  algorithm  is  used  to  estimate  a  Ham- 
merstein  model  based  on  simulation  data  of  a  SOFC.  The 
white-box  model  consists  of  a  set  of  coupled  equations,  as 
established  in  Section  2. 

The  inputs  to  the  system  are  the  temperature  input  (u\), 
the  fuel  flow  ( u2 ),  the  air  flow  (t/3),  and  the  current  (t/4).  The 
outputs  of  the  system  are  the  temperature  output  (y\)  and  the 
voltage  (y2). 

The  linear  subsystem  in  the  Hammer  stein  model  was  rep¬ 
resented  using  the  rational  orthonormal  bases  with  fixed  poles 
studied  in  Refs.  [17,22] 


Bi(q) 


{  q-!n  ) 


l  >  1 


(34) 


Bo(q)  = 


f  yi-itoi2\ 
V  q-%  0  ) 


(35) 


where  (§0,  §1,  •  •  %P- 1)  are  the  poles  of  the  bases. 

In  order  to  determine  the  model  order  of  the  linear  sub¬ 
system,  as  well  as  initial  guesses  for  the  location  of  the 
poles  of  the  bases,  the  same  input-output  data  were  used 
to  identify  a  linear  model  of  the  process  using  a  sub¬ 
space  method.  System  Identification  Toolbox  for  use  with 
Matlab  [23]  was  used  for  the  identification  of  the  linear 
model.  As  a  result  of  the  identification  process  a  fourth-order 
model  was  estimated  as  the  linear  part  of  the  Hammerstein 
model. 

On  the  other  hand,  a  third-order  polynomial  was  used 
to  represent  the  nonlinear  static  block  of  the  Hammerstein 
model.  This  choice  for  the  order  of  the  polynomial  was  a 
compromise  between  model  complexity  and  variance  error, 
and  the  decision  was  taken  based  on  the  values  of  Akaike’s 
Information  Theoretic  Criterion  (AIC). 

The  nominal  operating  conditions  of  the  SOFC  consid¬ 
ered  in  this  example  are  given  in  Table  1  [24,25].  For  the 


purposes  of  identification,  the  white-box  model  of  the  SOFC 
was  excited  with  band-limited  white  noise  around  the  nomi¬ 
nal  value  of  the  fuel  flow  ( u2 ),  while  all  the  other  inputs  were 
kept  constant  in  their  nominal  values.  Changes  in  the  fuel 
flow  were  produced  every  40  s,  with  a  maximum  amplitude 
of  ±100%  of  the  nominal  value.  A  set  of  10,000  data  was 
collected  from  the  simulation  with  a  sampling  time  of  0.1  s. 
The  first  5000  data  were  used  for  the  estimation  of  a  Ham¬ 
merstein  model  of  the  fuel  cell,  while  the  remaining  5000 
data  were  used  for  validation  purposes. 

The  true  (Full  model)  and  estimated  (Hammerstein  model) 
outputs  are  represented  in  Figs.  4-6.  A  step  change  in  the 
stack  current  (from  300  to  400  A)  is  applied.  The  true  (solid 
line)  and  estimated  (dashed  line)  output  voltages  correspond¬ 
ing  to  the  identified  model  are  represented  in  Fig.  4. 

Fig.  5  represents  the  cell  temperature  due  to  this  step 
change  in  the  stack  current. 

Finally,  the  voltage  step  response  of  the  SOFC  due  to  the 
change  in  fuel  flow  input  from  1 .2  to  1 .5  kmol  s-1  is  depicted 
in  Fig.  6. 


Fig.  5.  Output  stack  temperature  response  due  to  stack  current  step  increase. 
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Fig.  6.  Output  voltage  response  due  to  fuel  flow  step  change. 


5.  Conclusions 

This  paper  develops  a  multivariable  Hammerstein  model 
of  a  SOFC,  which  can  be  applied  for  small  signal  and  tran¬ 
sient  stability  studies.  The  model  based  on  electrochemical 
and  thermal  equations  describes  temperature  dynamics  and 
output  voltage  losses. 

Noniterative  algorithms  for  the  identification  of  this 
multivariable  Hammerstein  system  have  been  presented. 
The  key  issue  in  the  derivation  of  the  results  is  the 
representation  of  the  system  using  base  functions.  The 
suitability  of  the  proposed  methods  for  their  use  in 
the  identification  of  SOFCs  has  been  illustrated  through 
simulation. 


References 

[1]  J.  Padulles,  G.W.  Ault,  J.R.  McDonald,  An  integrated  SOFC  plant 
dynamic  model  for  power  systems  simulation,  J.  Power  Sources  86 
(1-2)  (2000)  495-500. 

[2]  DJ.  Hall,  R.G.  Colclaser,  Transient  modeling  and  simulation  of  tubu¬ 
lar  solid  oxide  fuel  cells,  IEEE  Trans.  Energy  Convers.  14  (3)  (1999) 
749-753. 

[3]  E.  Achenbach,  Three-dimensional  and  time-dependent  simulation 
of  a  planar  SOFC  stack,  J.  Power  Sources  49  (1-3)  (1994)  333- 
348. 

[4]  E.  Achenbach,  Response  of  a  solid  oxide  fuel  cell  to  load  change, 
J.  Power  Sources  57  (1-2)  (1995)  105-109. 

[5]  F.  Jurado,  Modeling  SOFC  plants  on  the  distribution  system  using 
identification  algorithms,  J.  Power  Sources  129  (2)  (2004)  205- 
215. 


[6]  K.  Narendra,  P.  Gallman,  An  iterative  method  for  the  identification  of 
nonlinear  systems  using  a  Hammerstein  model,  IEEE  Trans.  Autom. 
Control  11  (3)  (1966)  546-550. 

[7]  P.  Stoica,  On  the  convergence  of  an  iterative  algorithm  used  for 
Hammerstein  system  identification,  IEEE  Trans.  Autom.  Control  26 
(4)  (1981)  967-969. 

[8]  E.W.  Bai,  An  optimal  two-stage  identification  algorithm  for 
Hammerstein-Wiener  nonlinear  systems,  Automatica  34  (3)  (1998) 
333-338. 

[9]  M.  Boutayeb,  H.  Rafaralahy,  M.  Darouach,  A  robust  and  recursive 
identification  method  for  Hammerstein  model,  in:  Proceedings  of  the 
1996  IFAC  World  Congress,  1996,  pp.  447-452. 

[10]  F.  Chang,  R.  Luus,  A  non-iterative  method  for  identification 
using  Hammerstein  model,  IEEE  Trans.  Autom.  Control  16  (1971) 
464-468. 

[11]  T.  Hsia,  A  multi-stage  least  squares  method  for  identifying  Hammer¬ 
stein  model  nonlinear  systems,  in:  Proceedings  of  the  1976  CDC, 
1976,  pp.  934-938. 

[12]  S.  Billings,  S.  Fakhouri,  Identification  of  a  class  of  nonlinear  systems 
using  correlation  analysis,  Proc.  IEE  125  (7)  (1978)  691-697. 

[13]  W.  Greblicki,  Nonparametric  identification  of  Wiener  system,  IEEE 
Trans.  Inf.  Theory  38  (1992)  1487-1493. 

[14]  M.  Pawlak,  On  the  series  expansion  approach  to  the  identification  of 
Hammerstein  systems,  IEEE  Trans.  Autom.  Control  36  (6)  (1991) 
763-767. 

[15]  E.W.  Bai,  Identification  of  systems  with  hard  input  nonlinearities, 
in:  R.  Moheimani  (Ed.),  Perspectives  in  Control,  Springer,  Berlin, 
2001. 

[16]  R.  Haber,  H.  Unbehauen,  Structure  identification  of  nonlinear 
dynamic  systems — a  survey  on  input/output  approaches,  Automat¬ 
ica  26  (4)  (1990)  651-677. 

[17]  J.C.  Gomez,  E.  Baeyens,  Identification  of  block-oriented  nonlinear 
systems  using  orthonormal  bases,  J.  Process  Control  14  (6)  (2004) 
685-697. 

[18]  J.  Padulles,  G.W.  Ault,  J.R.  McDonald,  An  approach  to  the  dynamic 
modeling  of  fuel  cell  characteristics  for  distributed  generation  oper- 


152 


F.  Jurado  /  Journal  of  Power  Sources  154  (2006)  145-152 


ation,  in:  Proceedings  of  the  2000  IEEE/PES  Winter  Meeting,  vol. 
1,  2000,  pp.  134-138. 

[19]  M.D.  Lukas,  K.Y.  Lee,  H.  Ghezel-Ayagh,  Development  of  a  stack 
simulation  model  for  control  study  on  direct  reforming  molten  car¬ 
bonate  fuel  cell  power  plant,  IEEE  Trans.  Energy  Convers.  14  (4) 
(1999)  1651-1657. 

[20]  Fuel  Cell  Handbook,  fifth  ed.,  EG  and  G  Services  Parsons  Inc.,  U.S. 
Department  of  Energy,  2000. 

[21]  M.D.  Lukas,  K.Y.  Lee,  H.  Ghezel-Ayagh,  An  explicit  dynamic  model 
for  direct  reforming  carbonate  fuel  cell  stack,  IEEE  Trans.  Energy 
Convers.  16  (3)  (2001)  289-295. 


[22]  B.  Ninness,  F.  Gustafsson,  A  unifying  construction  of  orthonormal 
bases  for  system  identification,  IEEE  Trans.  Autom.  Control  42  (4) 
(1997)  515-521. 

[23]  L.  Ljung,  System  Identification  Toolbox,  Version  6.0.1,  The  Math- 
Works  Inc.,  Natick,  MA,  2004. 

[24]  J.A.  Kuipers,  in:  Proceedings  of  the  1998  Fuel  Cell  Seminar,  1998, 
p.  450. 

[25]  S.C.  Singhal,  Progress  in  tubular  solid  oxide  fuel  cell  technology, 
in:  Proceedings  of  the  1999  Solid  Oxide  Fuel  Cells,  Honolulu,  HI, 
USA,  1999,  pp.  39-51. 


